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Abstract 

Background: To determine tine optimal tliresliold of F-fluorodexyglucose F-FDG) positron emission 
tomograpliy CT (PET/CT) images tliat generates tine best volumetric match to internal gross target volume (IGTV) 
based on four-dimensional CT (4DCT) images. 

Methods: Twenty patients with non-small cell lung cancer (NSCLC) underwent enhanced three-dimensional CT (3DCT) 
scan followed by enhanced 4DCT scan of the thorax under normal free breathing with the administration of intravenous 
contrast agents. A total of 100 ml of ioversol was injected intravenously, 2 ml/s for SDCF and 1 ml/s for 4DCr. Then 

F-FDG PET/CT scan was performed based on the same positioning parameters (the same immobilization devices and 
identical position verified by laser localizer as well as skin marks). Gross target volumes (GTVs) of the primary tumor 
were contoured on the ten phases images of 4DCT to generate IGTViq. GTVpet were determined with eight different 
threshold using an auto-contouring function. The differences in the position, volume, concordance index (CI) and 
degree of inclusion (Dl) of the targets between GTVpfy and IGTViq were compared. 

Results: The images from seventeen patients were suitable for further analysis. Significant differences between the 
centric coordinate positions of GTVpEj (excluding GTVpEjisoJ and IGTViq were observed only in z axes (P < 0.05). 
GTVpETi5%, GTVpEr25% and GTVpet2.o were not statistically different from IGTViq (P < 0-05). G]Wp^^s% approximated closely 
to IGTViq with median percentage volume changes of 4.86%. The best CI was between IGTViq and GT\/p^j^s% (0.57). The 
best Dl of IGTViq in GTVpEj was IGTV,q in GTVpeti5% (0.80). 

Conclusion: None of the PET-based contours had both close spatial and volumetric approximation to the 4DCT IGTViq. 
At present 3D-PET/CT should not be used for IGTV generation. 

Keywords: Non-small cell lung cancer, Fluorodeoxyglucose positron emission tomography, Four-dimensional 
computed tomography. Standardized uptake value 



Background 

Worldwide, lung cancer is the most common cause of 
cancer-related mortality [1]. About 80% of the cases of 
lung cancer are non-small-cell lung cancer (NSCLC) [2]. 
Radiotherapy plays a major role in the management of 
patients with NSCLC who cannot tolerate or refuse 
surgery. Unfortunately, the prognosis of patients with 
NSCLC remains poor because of high rates of local 
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failure and distant metastases [3,4], with local control 
rates of approximately 50% after radical radiotherapy [5] . 
A geometric target miss induced by tumor motion 
during radiotherapy is considered as one of the main 
reasons for local failure [6]. In order to account for geo- 
metric uncertainties due to internal variations in tumor 
position, size, and shape, the International Commission 
on Radiation Units and Measurements (ICRU) report 62 
introduced the concept of an internal target volume 
(ITV) [7]. For lung tumors, respiratory motion is the 
major consideration for the ITV. 
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Currently, four-dimensional CT (4DCT) is widely used 
for the simulation of lung cancer. It is a reliable and 
effective tool for assessing tumor and organ motion [8,9] 
and can provide patient-specific information about 
tumor position, shape, and size at different phases of the 
respiratory cycle. The internal gross target volume 
(IGTV) can make the determination of the ITV more 
efficient [6]. IGTVio is generated by combining all 10 
individual gross tumor volumes (GTVs) contoured in 
each phase of the 4DCT dataset, which is thought to 
encompasses the motion information for the tumor in 
the whole respiratory cycle [10]. Therefore, it is expected 
to provide the most accurate IGTV based on a given 
4DCT dataset [11]. 

As a functional imaging modality, F-fluorodeoxyglu- 
cose (^^ F-FDG) positron emission tomography/computer 
tomography (PET/CT) images have been shown to have 
greater specificity and sensitivity than CT alone for the 
diagnosis and staging of NSCLC patients [12]. Further- 
more, the interobserver variability, as well as the intra- 
observer variability, could be significantly reduced when 
the F-FDG PET image was used for tumor volume de- 
lineation [13,14]. In addition, since three-dimensional PET 
(3D-PET) images are acquired over several minutes and 
represent the accumulated traces of multiple respiratory 
cycles, they may be capable of accounting for movement 
by indicating the average location of a tumor over time 
[15]. A phantom study by Caldwell et al. [16] concluded 
that PET imaging could more accurately depict the 3D 
volume of a moving phantom compared with spiral CT. 
Therefore, 3D-PET/CT might represent the ITV of a tu- 
mor. However, the optimal threshold values for patients 
with NSCLC have never been reported. 

Hence, We perform this study to investigate the 
appropriateness of the threshold method to determine 
the best volumetric match to 4DCT-based IGTVio when 
contouring the primary tumor volume of NSCLC. 
Additionally, the feasibility of 3D-PET/CT images was 
evaluated with respect to incorporating tumor motion 
into the radiation target for NSCLC. 

Methods 

Patients 

This study was approved by the Shandong Cancer Hos- 
pital and Institute review board, and 20 patients provided 
written informed consent. Patients with histologically 
proven primary NSCLC who were scheduled to undergo 
radiotherapy were eligible for this study, excluding those 
with atelectasis and/or obstructive pneumonia. None of 
them had previously been treated with radiotherapy or 
chemotherapy for their lung tumor. Between December 
2012 and December 2013, 20 patients were enrolled in this 
study. The maximal standardized uptake value (SUVmax) 
in the tumor of 3 patients was 2.89, 3.16 and 3.27, 



respectively. Therefore, they were not suitable for further 
analysis as several of the threshold-based contouring 
methods used would not discriminate between the tumor 
and background lung uptake. The other 17 patients in- 
cluded 13 men and 4 women, with a median age of 66 
years (range, 45-84 years). Six patients had centrally lo- 
cated lesions, and eleven patients had peripherally located 
lesions. The median of the SUVmax for the primary tumors 
was 11.34 (range from 6.07 to 25.51). Table 1 summarized 
the characteristics of the 17 patients and their primary 
tumors. 

CT simulation and image acquisition 

During the simulation, all patients were immobilized using 
thermoplastic mask for covering the head, neck and shoul- 
ders in the supine position. For each patient, an axial 
enhanced 3DCT scan of the thoracic region was performed 
followed by a enhanced 4DCT scan under uncoached free 
breathing conditions on a 16-slice CT scanner (Philips 
Brilliance Bores CT) with the administration of intravenous 
contrast agents. A total of 100 ml of ioversol was injected 
intravenously, 2 ml/s for 3DCT and 1 ml/s for 4DCT. De- 
tails of 3DCT and 4DCT scan as well as image acquisition 
were given in Li et al. [17]. Then, 3DCT and 4DCT images 
were transferred to MIM (MIM-6.0.4, MIM Software Inc, 
Cleveland, OH) imaging software. 

PET/CT simulation and image acquisition 

On the same day as the 4DCT scan, the FDG-PET/CT 
scans of the chest were performed with a integrated PET/ 
CT scanner (PhOips Gemini TF Big Bore). Through the 
same immobilization devices, the patient's position was 
identical to that for the 4DCT scan. Two radiation thera- 
pists were present to ensure the accuracy of the set-up by 
laser localizer and skin marks. All patients fasted for at least 
6 h before the PET/CT examination. All patients were 
injected with 7.4 MBq/kg body wei ght of ^**FDG and then 
rested for about 1 h in a quiet room before imaging. The 
16-slice CT component was operated with an X-ray tube 
voltage peak of 120 kV, 90 mA, a slice thickness of 5 mm 
and an interval of 4 mm, and was used both for attenuation 
correction of PET data and for localization of FDG uptake 
in PET images. No CT contrast agent was administered. 
PET scanning was performed covering the same axial range 
for 2 min per bed position (total of 3-5 bed positions). Both 
PET and CT acquisition was performed during free brea- 
thing. Data were reconstructed using an ordered subset ex- 
pectation maximization (OSEM) algorithm and attenuation 
correction derived from CT data. Then, the PET/CT 
images were transferred to MIM software. 

Image registration 

An initial automatic rigid registration was performed 
using MIM software. Due to the 3DCT and 4DCT 
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Table 1 Baseline characteristics of patients and the maximal standardized uptake values of every primary tumor 



Patients (n) 


Sex 


Age (y) 


Tumor location 


[Histology 


Tumor stage 




Tumor size* (mm) 


Tumor volume* (cm^) 


1 


M 


66 


Upper 


obe 


Adenocarcinoma 


T2N2 


25.51 


36 


17.49 


2 


F 


68 


Lower 


obe 


Adenosquamous carcinoma 


T2N2 


11 


46 


139.72 


3 


M 


70 


Upper 


obe 


Adenocarcinoma 


T2N2 


7.71 


18 


481 


4 


M 


79 


Upper 


obe 


Adenocarcinoma 


T1N2 


8.8 


33 


2640 


5 


F 


49 


Lower 


obe 


Squamous cell carcinoma 


T2N2 


8.85 


38 


38.59 


6 


M 


66 


Upper 


obe 


Squamous cell carcinoma 


T2N2 


1441 


33 


16.28 


7 


M 


75 


Upper 


obe 


Squamous cell carcinoma 


T2N2 


12.39 


32 


14.37 


8 


F 


65 


Upper 


obe 


Adenocarcinoma 


T2N2 


1341 


25 


8.12 


9 


M 


76 


Lower 


obe 


Squamous cell carcinoma 


T3N2 


15.13 


52 


1 1 542 


10 


M 


84 


Lower 


obe 


Adenocarcinoma 


T2N2 


24.81 


40 


27.53 


11 


M 


68 


Upper 


obe 


Squamous cell carcinoma 


T2N2 


9.83 


40 


45.97 


12 


M 


65 


Upper 


obe 


Adenocarcinoma 


T2N3 


14.38 


42 


39.43 


13 


F 


67 


Upper 


obe 


Adenocarcinoma 


T1N3 


6.99 


15 


5.96 


14 


M 


65 


Upper 


obe 


Squamous cell carcinoma 


T4N2 


14.52 


80 


157.67 


15 


M 


45 


Upper 


obe 


Adenocarcinoma 


TINS 


6.09 


10 


5.06 


16 


M 


65 


Lower 


obe 


Adenocarcinoma 


T1N2 


6.07 


24 


453 


17 


M 


60 


Upper 


obe 


Squamous cell carcinoma 


T4N1 


11.34 


34 


19.95 



SUVmax, maximal standardized uptake value; *Tumor sizes represent the major axis of the tumor; ^Tumor volume is the average volume of the 10 phases of the 
gross target volume delineated on four-dimensional CT. 



images for the same person were produced during the 
same imaging session, MIM would consider the images 
as being registered with each other. After the 3DCT and 
PET image datasets were co-registered with the help of 
the transmission CT from PET/CT, the 4DCT images 
would be auto-registered with the CT component of 
PET/CT. The registration was then manually adjusted 
by one radiotherapist, experienced in registering PET/ 
CT images, by matching bony anatomy such as the 
vertebral bodies at the level of the visible lung lesion. 
Hence, each contour was transferred to the 3DCT to 
calculate their specific parameters. 

Target volume delineation 

Our investigation focuses on the primary tumors. If the 
positive lymph nodes could not be separated from the 
primary tumor visually, they were delineated together as if 
they were part of the primary tumor. Patients were treated 
according to the 4DCT-based volumes and PET/CT con- 
tours were only used as part of a virtual planning study. 
Using the lung window setting (W = 1,600, C = -600) and 
mediastinal window settings (W = 400, C = 40) for the 
interface if the tumor was close to the mediastinum or 
chest wall [18], GTVs were manually contoured on all 10 
phases of the 4DCT scan by a single radiation oncologist 
and verified by another radiation oncologist. Both of them 
did not know the PET results in an effort to decrease bias. 
IGTVio were derived from the 10 phases of the GTVs. 
PET/CT-based GTV of the primary tumor (GTVpet) was 



defined by the auto-contouring function of MIM. After 
identification the primary tumor as a region of interest 
(ROI), MIM automatically calculated the SUVmax of ROI. 
Eight different threshold methods were used in this study: 
(1) SUV of 2.0 or greater (SUV2.0); (2) SUV of 2.5 or 
greater (SUV2.5); (3) 15% of SUV^^x within the ROI 
(SUV15%); (4) 20% of SUV„,ax within the ROI (SUV20%); 
(5) 25% of SUV^ax within the ROI (SUV25%); (6) 30% of 
SUVn^ax within the ROI (SUV30%); (7) 35% of SUV^ax 
within the ROI (SUV35%); (8) 40% of SUV„,ax within the 
ROI (SUV40%). All the noncancerous regions within the 
GTVpET) including the areas overlaid by the heart, bone 
and great vessels, were corrected to exclude manually with 
the help of the CT of component of PET/CT. 

Volumes comparison 

The differences in the position, size, concordance index 
(CI) and degree of inclusion (DI) between the GTVpEx 
and the IGTVio were compared. 

Target volume positions were defined by center of target 
coordinates and expressed using the x (left-right, LR), 
y (anterior-posterior, AP) and z (cranial-caudal, CC) 
coordinates of the center of mass. Centroid shifts in 
the 3D directions were calculated according to the 
formula as follows: 

3D centroid shifts = VlR^ + AP^ + CC^ 

The concordance index of volume A and B [CI (A, B)] 
was defined as the ratio of the intersection of A with B 
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variabilities of target volumes between GTVpet and 
IGTVio. We used the Spearman correlation test to 
analyze for associations between centroid shifts in the 
3D directions and CIs. Values of P < 0.05 were regarded 
as significant for all the tests. Descriptive statistics were 
used as appropriate. 

Results 

For lesion 15 and 16, the SUVmax in the tumor was 6.09 
and 6.07, we could not obtain GTVpeti5% target volumes 
as the volume obtained from the SUV15% contours were 
indistinguishable from background lung activity. 

Centroid shifts of the volumes derived from PET/CT and 
4DCT 

Table 2 showed the centroid shifts in 3D directions of 
the IGTVio volumes and the PET/CT volumes. The 
variations in the centroid coordinate positions in the 
CC direction of GTVpet2o% and IGTVio, GTVpet25% 
and IGTVio, GTVpet3o% and IGTVio, GTVpet35% and 
IGTVio, GTVpET40% and IGTVio, GTVpet2.o and IGTVio, 
GTVpET2.5 and IGTVio were significant (z = -2.131, 
-2.131, -2.012, -2.012, -2.012, -2.012, -2.226; P = 0.033, 
0.033, 0.044, 0.044, 0.044, 0.044, 0.026), while in the LR 
and AP directions were not significant (P > 0.05). The vari- 
ations in the LR, AP and CC directions of GTVpeti5% and 



Table 2 The centroid shifts of the GTVpet and IGTV,o in 3D directions (cm) 



Lesion 


GTVpETl 5% 


GTVpET20% 


GTVpET25% 


GTVpET30% 


GTVpET35% 


GTVpET40% 


GTVpET2.0 


GTVpET2. 


1 


1.10 


1.07 


1.04 


1.09 


1.08 


1.07 


1.08 


1.07 


2 


0.67 


0.75 


0.85 


0.94 


1.05 


1.14 


0.68 


0.79 


3 


0.77 


0.75 


077 


0.74 


0.70 


0.73 


078 


0.72 


4 


0.33 


0.34 


0.37 


0.38 


0.37 


0.35 


0.35 


0.38 


5 


0.38 


0.28 


0.23 


0.23 


0.25 


0.22 


0.28 


0.23 


6 


0.09 


0.16 


0.18 


0.18 


0.20 


0.20 


0.12 


0.14 


7 


1.13 


1.29 


1.36 


1.30 


1.43 


155 


1.16 


1.30 


8 


0.76 


0.75 


0.76 


0.75 


0.79 


0.77 


075 


0.73 


9 


0.22 


0.36 


0.48 


0.54 


0.59 


0.63 


0.14 


0.24 


10 


0.33 


0.29 


0.33 


0.32 


0.34 


0.35 


0.12 


0.23 


11 


0.20 


0.32 


0.36 


0.43 


0.46 


0.47 


0.31 


0.37 


12 


1.16 


1.14 


1.14 


1.11 


1.13 


1.17 


1.17 


1.14 


13 


0.32 


0.36 


0.41 


0.48 


0.49 


0.50 


0.44 


0.48 


14 


1.01 


1.02 


1.11 


1.17 


1.29 


152 


0.96 


0.98 


15 




0.15 


0.17 


0.25 


0.26 


0.22 


0.18 


0.42 


16 




2.01 


1.95 


2.02 


2.01 


2.01 


2.01 


1.98 


17 


0.43 


0.38 


0.43 


0.44 


0.50 


052 


0.33 


0.39 


Median 


0.43 


0.38 


0.48 


0.54 


0.59 


0.63 


0.44 


0.48 


sd 


0.37 


0.50 


0.49 


0.49 


0.50 


0.54 


0.51 


0.49 



sd, standard deviation. 



to the union of A and B [19]. The maximum value of CI 
is 1 if the two volumes are identical, and the minimum 
value is 0 if the volumes are completely nonoverlapping. 
That is, 

CI = ^ 
AuB 

The definition of DI of volume A in volume B [DI (A 
in B)] is the percentage of the overlap between volume 
A and B in volume A [20]. The formula is as follows: 

, AnB 
DI(AinB)=— . 

The DI can represent the percentage of one volume 
included by another volume, and 1-DI can represent the 
percentage of one volume not included by another 
volume. 

Statistical analysis 

Statistical analysis was performed using the SPSS 
software package (SPSS 17.0). The one-way ANOVA test 
was used to determine the variations in the DIs of 
GTVpET and IGTVio, and in the CIs of GTVpet and 
IGTVio- The Wilcoxon test was performed to estimate 
the differences of centroid coordinate positions between 
GTVpET and IGTVio, and also used to estimate the 
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IGTVio were not significant {z - 
P = 0.615, 0.378, 0.132). 



-0.502, -0.881, -1.505; 



Volume variation 

The volumes of primary tumors measured by PET/CT 
and 4DCT were shown in Table 3. Compared to IGTVio, 
GTVpETi5%; GTVpET2o% ^nd GTVpET2.o showed no signi- 
ficant difference {P values were 0.281, 0.102 and 0.687, re- 
spectively). Figure 1 illustrated the median percentage of 
volume changes from GTVpet to IGTVio standardized 
to the IGTVio for each case. The SUV 15% contour 
approximated most closely to the IGTVio with the lo- 
west median percentage volume changes of 4.86%. The 
corresponding values with respect to the IGTVio for 
the SUV20%, SUV25%, SUV30%, SUV35%, SUV40%, 
SUV2.0 and SUV2.5 contours were -8.87%, -24.51%, 
-37.51%, -43.57%, -52.24%, -7.03% and -21.93%, 
respectively. 

CI 

Figure 2 shows the median CIs of the PET volumes and 
the IGTVio volumes. The median CIs ranged from 0.30 to 
0.57 (F = 2.526, P = 0.018). The best fit for CI was at 
SUV15% (0.57), followed by SUV 2.0 (0.56). The CIs of 
IGTVio and GTVpet were inversely correlated with 
the centroid shifts in 3D directions (r = -0.668, -0.699, 
-0.728, -0.728, -0.801, -0.755, -0.711, -0.787; P = 0.007, 
0.002, 0.001, 0.001, 0.000, 0.000, 0.001, 0.000). 



Dl 

Figure 3 showed the median DIs of GTVpet in IGTVio, 
and IGTVio in GTVpet- The median DIs of IGTVio in 
GTVpET ranged from 0.31 to 0.80 (F = 7.814, P = 0.000). 
The best DI was IGTVio in GTVpeti5%- The median DIs 
of GTVpET in IGTVio ranged from 0.60 to 0.85 (F = 1.017, 
P = 0.422). The best DIs were GTVpet35% and GTVpet4o% 
in IGTVio. 

Discussion 

Accurate definition of the target volume, ideally incor- 
porating metabolic information, becomes paramount 
importance in the current trend in NSCLC treatment 
planning [21]. A number of studies compared 3DCT vol- 
umes with ^**F-FDG PET/CT volumes of NSCLC [22-24]. 
However, the best methodology for applying F-FDG 
PET/CT to IGTV definition is not currently well estab- 
lished. To the best of our knowledge, there is few study 
that has compared tumor sizes and CI values between 
GTVpET and IGTVio in contouring NSCLC. 

In this study, we determined the IGTVio from 10 phases 
of the 4DCT dataset and used them as the reference to 
find the optimal threshold that yield the best match bet- 
ween the GTVpET and IGTVio in both the target size and 
the spatial conformity. Our study revealed that GTVpet 
using a threshold setting of SUV15% approximated most 
closely to the IGTVio with the lowest median percentage 
volume changes. When using the threshold level of 



Table 3 The volumes of primary tumors measured by PET/CT and 4DCT (cm^) 



Lesion 


GTVpET15% 


GTVpET20% 


GTVpET25% 


GTVpET30% 


GTVpET35% 


GTVpET40% 


GTVpET2.0 


GTVpET2,5 


IGTVio 


1 


24.98 


20.59 


17 


15.31 


13.83 


11.94 


32.75 


30.37 


24.5 


2 


159.69 


122.17 


95.89 


72.35 


54.9 


40.9 


135.47 


107.23 


173.72 


3 


16.34 


10.47 


7.82 


6.13 


4.15 


3.59 


7.67 


5.04 


5.9 


4 


39.49 


34.32 


29.54 


26.89 


23.72 


21.54 


31.11 


28 


37.66 


5 


70.9 


53.99 


46.02 


41.27 


3648 


32.39 


4945 


4245 


44.01 


6 


19.41 


14.44 


114 


9.04 


8.07 


7.28 


2046 


16.81 


22.31 


7 


19.68 


13 


8,5 


5.93 


4.27 


3.52 


18.07 


12.76 


19.99 


8 


13.03 


10.16 


7.9 


649 


5.53 


4.68 


13.03 


10.70 


9.8 


9 


183.15 


140.37 


115.35 


98.64 


8646 


76.09 


201.56 


1 66.58 


152.8 


10 


24.24 


19.76 


1635 


13.81 


11.99 


10.33 


39.27 


33.15 


33.16 


11 


72.69 


56.06 


46.42 


39.95 


3447 


29.97 


54.44 


45.79 


56.93 


12 


45.43 


35.7 


28.19 


2345 


18.3 


15.83 


4846 


4040 


57.53 


13 


26.49 


19.64 


14.06 


9.92 


748 


57 


10.84 


7.15 


117 


14 


141.98 


112.7 


89.36 


69.05 


51.05 


35.97 


152.57 


128.12 


164.1 


15 




8.59 


5.76 


3.94 


3.12 


2.05 


5.91 


3.67 


7 


16 




7.64 


5.68 


4,2 


3.55 


2.84 


3.78 


2.66 


11.82 


17 


42.34 


31.21 


27.05 


23.24 


19.18 


16.89 


33.02 


28.60 


26.15 


Median 


39.49 


20.59 


17 


15.31 


13.83 


11.94 


32.75 


28.60 


26.15 


P value 


0.281 


0.102 


0.004 


0.000 


0.000 


0.000 


0.687 


0.031 





P value, GTVpET VS IGTV,o, 
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Figure 1 Median percentage volume change of GTVpEj to 
IGTVio. SUV n = SUV of n; SUV n% = n% of maximum SUV 



>SUV25% and/or > SUV2.5, the PET-based tumor sizes 
were estimated to be smaller than the IGTVio. Therefore, 
on the basis of the results of our study, the SUV threshold 
setting of >25% and/or >2.5 is not suitable for IGTV 
contouring in NSCLC. 

Analogously, Hanna et al. [25] compared volumes 
from a manual method and five automated PET segmen- 
tation techniques to 4DCT-derived ITV and found that 
none of the PET target volumes approximated closely to 
the 4DCT target volumes. However, in their study, the 
patient's PET/CT and 4DCT scans were not acquired on 
the same day or in identical position. In this circum- 
stance, it was possible that changes in tumor geometry 
or size occurred and potentially increased the likelihood 
of a mismatch between the PET-based contours and the 
4DCT-based contours [25]. Caldwell et al. [16] reported 
that using a threshold as low as 15% of the maximum 
value could account for respiratory motion and more ac- 
curately depict the true extension of the moving target. 
Another phantom study by Okubo et al. [26] concluded 
that when a threshold value of 35% of the measured 
maximum FDG activity was adopted, the sizes of PET 
delineation were almost the same for static and moving 
phantom spheres of 22 mm or more in the axial plane. 
Our study was similar to the result of Caldwell et al., but 
smaller than the result of Okubo et al. This is possible 
because patients enrolled in our study had a range of 
tumor sizes and positions. Moreover, the 15% threshold 
method was not suitable for contouring some lung 



tumors that have low SUV, because it might fail to dis- 
tinguish tumor from background lung activity. Neverthe- 
less, Okubo et al. suggested that the threshold of 35% of 
measured maximum FDG activity was only a provisional 
criterion for tumors of 2-4 cm given that appropriate 
threshold values could be changed on the basis of the 
tumor size [26]. In addition, it should be acknowledged 
that any phantom studies versus clinical comparison is 
limited. For example, unlike in real tumors, the FDG 
distribution in the spheres of the phantoms and in the 
background was homogenous. 

Similarity in absolute volume does not mean identity in 
the space location. Our results indicated that GTVpeti5%) 
GTVpET2o% and GTVpet2.o showed no significant dif- 
ference with IGTVio in target volume. However, the CIs 
of them were significantly lower than 1.0. The best CI was 
between IGTVio and GTVpeti5%> which was only 0.57. It 
is not surprising as GTVpeti5% is the biggest volume and 
hence has the greatest degree of potential overlap. Based 
on this consideration, this may not make it the most 
accurate. The poor CIs suggested great unconformity bet- 
ween what was indicated abnormal on PET image and on 
CT image. One of the reasons is that shape and/or pos- 
itional alterations between IGTVio and GTVpet had oc- 
curred. Our study showed that the minimum variation 
in the centroid coordinate position in the 3D direction 
of GTVpET and IGTVio was 0.38 cm (median). More- 
over, the CIs were inversely correlated with the centroid 
shifts in 3D directions. Although patients in our study 
were immobilized in the same position for both 4DCT 
and PET/CT, millimetric precision in set-up using 
immobilization devices may not be feasible. Further- 
more, a rigid registration might not be sufficient for 
lung tumors. Hence, registration error may inevitable 
affect the spatial position between GTVpet and IGTVio- 
In addition, it is possible that some of this difference 
may be related to differences in the patient's breathing 
pattern between acquiring the PET/CT and 4DCT. Dif- 
ferent breathing pattern can influence tumor size, shape 
and distribution of activity on the free-breathing PET 
images [27]. 
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Hanna et al. [25] used the Dice similarity coefficient 
(DSC) to assess volumetric, shape and positional simi- 
larity in the PET-generated target volumes and 4DCT 
target volumes. Their study revealed that the highest 
DSCs (mean) were 0.64. Grills et al. [21] investigated 
the impact of PET/CT for GTV definition in NSCLC 
using the matching index (similar to CI in our study) 
to compare the GTV, as defined by CT, with the GTV 
defined by PET co-registered with CT. In their study, 
the mean matching index was 0.65. Gondi et al. [23] 
demonstrated that CI values of NSCLC with the 
incorporation of FDG-PET and CT were 0.44. The 
results of our study were similar to data published in 
prior studies [21,23,25]. In addition, although CI can 
provide the most information on both volume change 
and positional change [28], it cannot quantify the 
percentage of one volume included by another volume. 
Further analyzing the inclusion relation between 
GTVpETi5% and IGTVio, there was 20% of (median) 
GTVpETi5% not included in IGTVio and 40% of (median) 
IGTVio not included in GTVpeti5%- It suggested that 
ITGVio did not encompass GTVpeti5% completely or vice 
versa. Therefore, we concurred with Gondi [21] who 
concluded that although the quantitative absolute target 
volume could sometimes be similar between CT and 
PET, the qualitative target locations can be significantly 
different. 

One limitation of our study is the small number of 
patients studied, so that subgroup analyses were limited. 
Therefore, a larger cohort of patients with many differ- 
ent tumor sizes and locations should be conducted to 
further investigate the relationship of using 3D-PET/CT 
and 4DCT for contouring IGTV of NSCLC. We are 
continuing our work to enroll more patients for further 
clinical investigations. 

Conclusion 

None of the PET based contours had both close spatial 
and volumetric approximation to the 4DCT IGTVio. At 
present 3D-PET/CT should not be used for IGTV 
generation. 
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